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Abstract 

We investigate the infinite volume limit of quantized photon fields in multi- 
mode coherent states. We show that for states containing a continuum of coherent 
modes, it is mathematically and physically natural to consider their phases to be 
random and identically distributed. The infinite volume states give rise to Hilbert 
space representations of the canonical commutation relations which we construct 
concretely. In the case of random phases, the representations are random as well 
and can be expressed with the help of ltd stochastic integrals. 

We analyze the dynamics of the infinite state alone and the open system dy¬ 
namics of small systems coupled to it. We show that under the free field dynamics, 
initial phase distributions are driven to the uniform distribution. We demonstrate 
that coherences in small quantum systems, interacting with the infinite coherent 
state, exhibit Gaussian time decay. The decoherence is qualitatively faster than 
the one caused by infinite thermal states, which is known to be exponentially rapid 
only. This emphasizes the classical character of coherent states. 


1 Introduction 

Coherent states play an important role in the theory and experiment of quantum me¬ 
chanics. They have first been discovered by Schrodinger in 1926 as quantum states 
which behave in many respects like classical ones [17]. Fourty years after, Glauber [7,8] 
realized that these states are particularly suited to describe optical coherence, hence 
their name. Their use in modern quantum optics is now ubiquitous [6,11,18] and more 
recently, coherent states have been used successfully in quantum information experi¬ 
ments (for instance in implementations of quantum key distribution [9]). 

In this paper, we analyze the infinite volume limit, or thermodynamic limit, of 
the quantized radiation field in a multimode coherent state. We start with the field 
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confined to a box A C The field modes (or, momenta of the associated particles) 
are quantized and take only discrete values. To each discrete mode k is associated 
a creation and an annihilation operator a|. and a^, respectively, satisfying the usual 
canonical commutation relations [afc,a|] = Sk/ (Kronecker symbol). To each mode k 
corresponds as well a family of coherent states \a)k, indexed by a G C. They are defined 
to be normalized eigenvectors of the annihilation operator o^, satisfying ak\a)k = Oi\a)k- 
The state | a) k does not have a definite number of particles as it is not an eigenvector of 
the number operator N = The average of N in \a)k is |ap. We set a = \a\e'^ 

and call 9 the phase of the coherent state. The state \a)k can be obtained by applying 
the displacement operator Dk{a) = to the vacuum vector n of the quantum 

field, \a)k = Dk{a)Q. This is a single mode coherent state. If one selects N (discrete) 
modes ki,..., k^, the associated multimode coherent states are ^^^(ai) • • • Dkj^{aj\f)^, 
for any choice of G C, j = 1,..., A^. 

It is well known that the electromagnetic radiation generated by a classical current 
is a multimode coherent state of this form, and so is the light produced by a laser in 
certain regimes [11,12]. One of the physical motivations to consider the infinite vol¬ 
ume limit is that the radiation quantum field provides an environment (“reservoir”) 
which is extremely large, spatially, relative to objects of interest that are placed in 
this environment. “Object-environment systems” are called open quantum systems. 
The archetypical example is a spin 1/2 (a “qubit”) interacting with an environment 
(photons, phonons). The infinite nature of the reservoir causes irreversible dynamical 
effects in the spin, such as thermalization, decoherence and disentanglemet (in pres¬ 
ence of several spins). On the mathematical side, taking the infinite volume limit is 
interesting in its own right, as it uncovers new Hilbert space representations of the 
canonical commutation relations. Indeed, in 1963, Araki and Woods [1] considered 
the infinite volume limit of “Fock states” having the form a*{fi) ■ ■ ■ a*{f]\f)Q, which 
describe thermal equilibrium as well as condensate states of the field. They found the 
famous Araki-Woods representations of the canonical commutation relations. These 
representations, as well as their fermionic counterparts, the Araki-Wyss representa¬ 
tions [2], have proven to be an important tool in the mathematical analysis of open 
quantum systems close to thermal equilibrium. 

In an effort to be able to handle reservoirs, that is, spatially infinitely extended 
systems with non-vanishing particle density, which are physically different from thermal 
ones, we tackle here the question of coherent reservoirs. As mentioned above, they 
describe the electromagnetic field created by classical sources or by lasers. 

Outline of main results. We start off with the quantized radiation field in a 
finite box A C M^, pick modes ki,... ,k]\f and consider the multimode coherent state 
Dki (cq) • • • Dkj^ (a 7 v)fl- What happens to this state when the size of the box increases, 
A —)■ ? As A changes, so do the eigenmodes of the quantum field, and in the limit of 

infinite volume, the values of the modes becomes a continuum, /c G Any state u of 
the radiation field (in finite or infinite volume) is determined uniquely by its expectation 
functional E{f) = uj{W{f)), where W{f) is the Weyl operator smoothed out with a 
test function /. Convergence of a sequence of states (with A„ —)■ R^) is then meant 
as convergence of the associated En{f), for all /. 
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(A) Thermodynamic limit. We look at two basic scenarios. 

• Infinite volume limit for N fixed modes: Fix ki,... ,k]\f G and consider the 

limit A —)• while keeping the particle densities pj = |q;jP/|A|, and the phases 

9j ^ , j = 1 , ..., N , fixed. 

We show that the A-mode state (i.e., the associated expectation functional) has a 
thermodynamic limit and calculate the latter explicitly in Proposition 2.1. The result 
depends, of course, on the densities pj and the phases 6j, j = 1,... N. 

• Infinite volume limit and continuous mode limit: Fix a mode density distribution 
p{k), meaning that p{k)d^k is the spatial density of particles (number of particles 
per unit volume in direct space) having momenta in an infinitesimal volume d^k 
around k. Consider the infinite volume limit of the state associated to p{k). 

The “naive” infinite volume and continuous mode limit of the multimode coherent state 
(i.e., its expectation functional) does not exist, see (2.10). However, the state obtained 
by mixing the phases of the N modes in an identical and independent way according 
to a measure p on the unit circle, does converge. The limit depends on the mode 
density distribution p{k) as well as on the phase distribution p, see Proposition 2.2. It 
is interesting to note that in laser experiments, the intensity of the laser field (measured 
by p{k)) is observed to be almost free of fluctuations due to saturation properties of 
the laser. However, the phases are subject to fluctuations and drift randomly. The 
theoretical framework describing these properties is called the randomly phased laser 
model (see e.g. [11], paragraph 11.8.6). The fluctuation of the phases is precisely what 
makes the infinite volume and continuous mode limit exist in our analysis. Considering 
the phase mixed state is the same as taking an expectation of random phases, and 
the above mentioned result (Proposition 2.2) is equivalent to the convergence of the 
expectation of the state with random phases. We then ask if a stronger convergence 
holds. Consider the phases to be iid, distributed according to a measure p on 5^. 
We show in Proposition 2.3 that the corresponding random phase state converges in 
distribution. The convergence is in essence due to the central limit theorem. The 
random part of the limit expectation functional E{f) is given by where xif) 

is a complex valued random variable which we construct as a suitable Ito stochastic 
integral. It has the property Rex(/) ~ W(0, cj^(/)^), where the latter is a (real) 
normal random variable with mean zero and variance (T^(/)^ which depends on the 
phase distribution p. This result is the content of Theorem 2.4. We point out that 
only the second Fourier moment p{2) (c.f. (2.16)) of the phase distribution p enters 
the infinite volume expression. 

(B) Hilbert space representation. In Section 2.4 we construct explicit Gelfand- 
Naimark-Segal (GNS) Hilbert spaces of the three infinite volume states constructed be¬ 
fore, namely, the finite modes state, the phase averaged state and the random phases 
state. Those representations are given in Theorem 2.5. They all are regular repre¬ 
sentations of the canonical commutation relations, hence defining (represented) field-, 
creation- and annihilation operators (Proposition 2.6). 
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(C) Dynamics of the infinite coherent state. The Heisenberg dynamics of the 

field is a Bogoliubov transformation W{f) i—)■ where t is time and (in Fourier 

space) e = £{k) is the dispersion relation. For the infinite volume state with N fixed 
coherent modes kj and phases 6j, the effect of the dynamics is to rotate the phases, 
6j I—)■ 6j —te{kj), making the infinite N mode coherent state a quasiperiodic function of 
time. However, the situation of the mixed phase state and of the random phases state is 
quite different. We show in Proposition 2.8 that the dynamics drives those states into a 
final state, which corresponds to the uniform distribution of the phases. More precisely, 
starting off with the infinite volume state with phase distribution ^ on the circle (and 
continuous mode distribution p), as time moves on, the system converges to the infinite 
volume state with phase distribution dp(0) = d0/27r (and the same continuous mode 
distribution p). In this sense, the uniform phase distribution is the stable one. 

(D) Coupling to an open quantum system. In Section 2.5.2 we consider an 
A^-level quantum system in contact with the infinite volume coherent reservoir having 
uniformly randomly distributed phases (the dynamically stable reservoir state as ex¬ 
plained in point (C) above). The coupling between the Wlevel system and the reservoir 
is energy conserving, meaning that the interaction term in the Hamiltonian describing 
the coupled evolution commutes with the system Hamiltonian. Such couplings are often 
considered to describe so-called phase decoherence without energy exchange [10,13,16]. 
This model is explicitly solvable in the sense that we can calculate the exact density 
matrix of the Wlevel system at all times (c.f. (2.55), (2.57)). We show that the expec¬ 
tation of the off-diagonal density matrix elements (in the energy basis) have Gaussian 
time decay due to the coupling with the coherent reservoir. This is a striking difference 
relative to the thermal reservoir case, where this decay is only exponential. We conclude 
that a small system placed in a coherent reservoir undergoes much faster decoherence 
than in a thermal environment.^ This very rapid loss of “quantumness” (encoded by 
coherence of the small system) is yet another manifestation that coherent states behave 
to some extent classically. 


2 Setup and main results 

We consider non-interacting quantum particles confined to a box of sidelength L in d 
dimensions, 

A = [-^2,^2]^^ C 

The wave function of a single particle is an element of T^(A, dx) with periodic boundary 
conditions. The space of pure states of the system of particles is the symmetric Fock 
space [4] 

F^F{L\k,dx))=^L%^^{K-,d^x). 

n>0 

^The name “coherent states” refers to the quantum field and is motivated by the fact that correlation 
functions of the field factorize in those states, which is the same as for classical coherent fields [11]. 
On the other hand “decoherence” of a quantum system is an entirely different notion, which refers to 
a system losing quantum correlations and becoming close to a classical one [10]. 
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A state is then given hy ip = ©n>oV'^”'\ where • • •, Xn) is a symmetric function 

of n variables Xj G with periodic boundary conditions in A”. The summand n = 0 
of the Fock space is called the vacuum sector, it equals C and is spanned by the vacuum 
vector n (such that 17 = 1 , = 0 for n > 1 ). 

The dynamics obeys the Schrodinger equation 

'ipt = 

where the self-adjoint Hamiltonian is the second quantization of a one-body Hamilton 
operator. For photons (massless relativistic particles), H is the square root of the 
Laplace operator, 

n 

{Hlp)^'^\xi, . . . ,Xn) = ^ (Xi, . . .,Xn), 

i=i 

understood as a self-adjoint operator with periodic boundary conditions. 

The creation operator a*{f) is defined for / G T^(A, dx) by 

..., Xn) = VnSf{xi)il>^'^~^^ (X 2 , • • •, Xn), 

where S is the operator of symmetrization over the variables xi,... ,x„. The adjoint 
of the creation operator is the annihilation operator a{f), given by 

(a(/)V’)^”^(xi,... ,Xn) = \/n"+T f /(x)V'("-+^^(x,xi,... ,Xn)dx. 

Ja 

The self-adjoint field operator is defined by 

and the unitary Weyl operators is 

TT(/) = 

It is practical to introduce the operator valued distributions a*(x), a(x) by setting 

a*if)= [ fix)a*{x)dx, a{f)= f /(x)a(x)dx. 

Ja Ja 

The canonical commutation relations take the following equivalent forms: 

[a{f),a*{g)] = {f,g) 

[a{x),a*{y)] = d{x - y) 

[$(/), ^>( 5 )] = ilm {f,g) 

W{f)W{g) = e-ii“<^’^>IT(/ + 5 ). (2.1) 

Let V' be a normalized vector in iF. It determines an expectation functional, defined 
by E{f) := {ip,W{f)ip), / G L^(A,dx). Conversely, any functional E : L^(A, dx) —)• C 
satisfying the three conditions 
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(El) E{0) = 1 
(E2) 'm) = E{-f ) 

(E3) Ylk,k'=i - /fc/) > 0, for all K > 1, Zfc G C, G L2(A,dx) 

determines a state uj on the C*-algebra generated by the Weyl operators by the relation 
u;{W{f)) = E{f), see for instance [14], 

To define coherent states of the particles in the volume A it is convenient to pass 
to the momentum space representation, the Eourier transformation of the Eock space 
E. There coherent states take a simple form, while their expression in position space 
(as vectors in E) is more cumbersome. 


2.1 Momentum space representation 


The single-particle Hilbert space L^(A, dx) is unitarily equivalent to via the 

Eourier transform ^ : L^(A, dx) —)■ 


idfm = fk = T-"/' / e-'^"/(x)dx, (2.2) 

Ja 

having the inverse 

Here, kx is the dot product k ■ x and the factors guarantee that ^ is unitary. 

Accordingly, the Fock space E is unitarily equivalent to its momentum version 

jt ^ E{l\^Z'^)) = ©„>o (/'(f 

Let H and fl, be the vacua of the Fock spaces E and E, respectively. The unitary map 
between the Fock spaces is given by the natural lifting of 5^, 

• • • a*{fi)n = a*ih) • • • a*(/^)H. 

Accordingly, the creation operators transform as = a*{f). We write 

aif)= E 

k&'^'&d- 


and so 


«•(/) 





dx. 


Comparing this with 5^0* (/)5^ ^ 


/(x)5^a*(x)(J ^dx yields the relations 


k&2^Zd 


5-I4 (5 = [ e''^^a*{x)dx. 

Ja 


The field- and Weyl operators are transported to the momentum space as 

= Hf) and = w{f), 

where 4'(/) = ^ + hak) and W{f) = 
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2.2 A^-mode coherent states in finite volume 

The coherent state associated to the collection of N modes k[,... ,k'j^ G and N 

complex numbers ai,..., otv is the normalized vector 


T = e 


E 




(2.3) 


The expectation functional of the coherent state (2.3) is defined by 

= (2.4) 

for all / G The operator a'^, ay, is the number operator of the mode k'^ G 

Its average in the coherent state is 



and can be interpreted as the intensity of the mode in question. 


2.3 Infinite volume and continuous mode limits, random phases 
2.3.1 Infinite volume 

The momenta in the finite-volume coherent state are given by /c'- = 27rnj/L, for nj G Z'^. 
As L increases, the spacing of the momenta becomes increasingly small. Let now 
/ci,..., /cat G be A" arbitrary (‘continuous’) momenta and let nj = nj{L) G Z'^ be 
such that kj{L) = 27rnj{L)/L satisfies limA^oo kji^) = kj, j = 1,..., N. We want to 
take the thermodynamic limit of (2.4), 

hm E^if) = Eyif). (2.5) 

L^oo 

This means we take /c' = kl{L) and L —)■ oo, while keeping fixed the particle densities 
Pj > 0 which count the (average) number of particles in mode k'j per unit volume, for 
j = 1, ..., A. In other words, we have |ajP = L^Pj, or 

a,(L) = L‘'/2^e'"L (2.6) 

where 9j is the phase of the complex number aj. 

Proposition 2.1 (Thermodynamic limit for A modes) Let ki,...,k]sf G M and 

Pi,..., Pn > 0 be arbitrary momenta and arbitrary particle densities and suppose that 
f G L^(M‘’*,dx) n L^(M‘’*,dx). Then the limit (2.5) exists and 

EnU) = i?Fock(/) (2.7) 

where Epockif) = e"2ll^ll^ and f{k) = e"'*’"’/(x)dx. 
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Here, -EFock(/) is the Fock expectation functional, determined by the vacuum state, 


EfoM = = e--M = e-2(2-)-‘*ll/li. 


Here, || • II 2 is the L^-norm (of functions of A: G or x G W^), 

II/II2 = [ l/(a;)|^dx and ll/lli = f \f{k)\^dk. 

JRrf JRd 

Remark. By adopting the definition f{k) = e“‘^*/(x)dx of the Fourier trans¬ 
form (see Proposition 2.1), we obtain /(x) = (27r)“‘^ e'^^/(A:)dA: and 




2.3.2 Continuous modes 

We are now interested in an infinite-volume coherent state which contains a contiuum 
of modes. One may perform the infinite-volume limit and the continuous mode limit 
simultaneously, or one can take the continuous mode limit of (2.7). The result is the 
same and we do the latter (see Section 3.4). Let p{k) be a prescribed mode density 
distribution. That is, given a cube I C M'^, the integral fj p{k)dk is the spatial density of 
particles in the infinite volume state, having momenta in I. Consider p to be supported 
in a finite cube [—ii, Rff. Discretize the cube by taking an N (large) and setting 

kj = (-i? + ..., -i? + j^^) G (2.8) 

where ji,..., ^ {1,2,..., A^}. Here, we view j = {ji,... ,jd) as a multi-index. 

In the previous considerations, pj was the number of particles per unit volume hav¬ 
ing momentum kj. Since p{k) is, by definition, the number of particles per unit volume 
in space and per momentum volume dk, we have p{kj) = with Akj = {2R/N)'^. 
Consequently, the sum in the phase of the infinite volume expectation functional (2.7) 
is 


e-'^<^^^J{k,) = {2R/Nf^ Y, e-'^^^^^^,j2p{k,)f{k,). (2.9) 

Here, 6{h) is a function we can choose, which determines the phase of the mode k. For 
N large, (2.9) equals approximately 

{N/2Rf/^ [ f{k)dk ~ (2.10) 

which diverges as N ^ 00 . It follows that the infinite volume discrete mode expectation 
functional EN{f), (2.7), does not have a continuous mode limit in this simple sense. 



2.3.3 Phase mixture 


Suppose the phases in the finite-volume coherent state (2.3) are not fixed, but dis¬ 
tributed independently according to probability measures dfij on the circle. For given 
phases, denote the coherent state by The mixed state is given by the density 

matrix 


r27r 


C27T 


Q = 


d^i(0i) 


dllN{0N) 


It follows immediately from Proposition 2.1 that the associated expectation functional 
(/) = Tr(^IF(/)) has the infinite volume limit 


^ /‘2'k 


{En) if) = £;Fock(/) n / (2.11) 


Here, we use the notation ( ) to indicate that we have taken the average over phases 
(phase mixture). For identically and uniformly distributed phases, d/xj = (27r)“^d0, 
(2.11) becomes 


N 


(En) if) = i?Fock(/) n \fikj)\), 


where 


i=i 

27r 


Jo(Va2 + 62) = f " ^Q-Kaco,e+b,me) 

Jo 27 r 


( 2 . 12 ) 


(2.13) 


is the Bessel function. 

The continuous mode limit of the uniformly distributed, phase averaged expectation 
functional is now well defined. To see why this is the case, consider the phases 
6j = 0{kj), where kj is given in (2.8). The averaged functional is (see also (2.9)) 

i(2iJ/Ar)'i/2Ree-iV2p(fci) f(kj) 


{En) if) = i^Fock(/) n 


27r 


Since 


p27r 


de 


—e“‘^ = 0 
27r 


(2.14) 


the first order term in the expansion of the integral in powers of N vanishes, and 
we have 


{E^) if) = F;Fock(/) n 

je{i,...,N}<^ L 


j _ (2ityp(k,)\f{k,f ^ 


(2.15) 


Consequently, 

N^oo Epockil) 


-1 hm i2R/N)‘^ Y. + OiN-^/^) 


= -h j P(.k)\fik)\^dk. 
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This derivation, which we have carried out for / and p continuous and compactly 
supported, can be extended to / G by a density argument. 

In case the phases are identically and independently distributed according to a 
measure d^ on the circle, the continuous mode limit exists if and only if p{l) = 0, 
where 

r27r 

pin) := / d^(0)e“'"'^, n G Z. (2-16) 

JO 

This is so because ^(1) = 0, which is the analogue of (2.14) above, is equivalent to the 
vanishing of the term in (2.15). A little generalization of the above calculation, 

where d^ was uniform, shows the following result. 

Proposition 2.2 Let p{k) he a eontinuous distribution of momenta per unit spatial 
volume, having compact support. Suppose the phases are identically and independently 
distributed on [0,27r], aecording to a probability measure p satisfying p{l) = 0 (see 
(2.16)J. Then the phase mixed expectation functional (2.11) has the continuous mode 


limit 


hm (Sw) if) = iE) if) = ^Fock(/)e-2'^-(^) , 

N^oo 

(2.17) 

where 

/* 


^uiff ■= y^^/5(^)(l/(^)P + Re{/l(2)/(A:)2}jd/c. 

(2.18) 

2.3.4 

Random phases 



Considering the phases 9j = 6j{io), j G {1,... , to be random variables, the 
coherent state T = given in (2.3) is a random pure state defining the random 
expectation functional (in infinite volume, by Proposition 2.1) 

En,M) = (2.19) 

where 

e,(w) = f{k,). (2.20) 

Denote by IE the expectation over the randomness (tc). If the phases are independent 
and identically distributed over [0, 27r] according to p with ju(l) = 0, then Proposition 
2.2 states that 

hm E[En,M)] = (E) if). 

N^oo 

On the other hand, we know from (2.10) that iIlAr,a;(/) does not converge almost every¬ 
where. Does E]\}^^if) converge in a sense that lies in between these two? The answer 
is given by the central limit theorem. 

Proposition 2.3 Suppose the phase distribution satsifies /i(l) = 0. Then 

E ei(cj) A A?^(0,a^(/)2), asN^oo, (2.21) 

The right hand side is a normal random variable with mean zero and varianee (j‘fi{f) 
given in (2.18). The convergence is in distribution. 
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The meaning of convergence in distribution is that the distribution function 


F^(x)=P ^ 

\ je{i,...,N}d 

converges pointwise to that of the normal with mean zero and variance cr^{f)‘^, 
lim F]\f{x) = —=^ - [ 

for all X G M. 

It follows from the fact that x e'* is bounded and continuous and (2.19) that 


En,M) ^ £;,.(/) = asAT^oo. (2.22) 

Our task now is to hnd an explicit representation of the random variable A/’aj(0, it^(/)^) 
such that the corresponding (2.22), defines an expectation functional. Let Xi^{f) 

be such a random variable. It is easily seen that E^^{-) satisfies conditions (E1)-(E3) 
provided that X^{-f) = -X^{f) and X^{fi + / 2 ) = X^{fi) + X^{f 2 ). We are then 
taking X^{f) = ReXcoif), where, for / G L^(M^,dx), we define Xlo as an Ito integral 


with 


xM)= [ dBf{k)Si{k)m + l [ dB^{k)S 2 {k)m, (2.23) 

jRd Jm<i 


Siik) = 

1 P(.k) 

(l + (i(2)) 

(2.24) 

Y 1 + Re/i(2) 

S2{k) = 

/ P(.k) 

yi - iP(2)p. 

(2.25) 

Y f T Re/i(2) 


Here, Bf and B 2 are two independent Brownian motions of dimension d, as we are 
going to explain below. 

Remarks. The choices of Si, S 2 are not unique. Since |/i(2)| < 1 one easily sees 
that both S'!(A:) and S 2 {k) are bounded as functions of /i(2). If Re/i(2) = — 1, then one 
may take Si = iy^ and S 2 = y/p- This happens e.g. if dp, = \{dT ^/2 + <^- 77 / 2 }- 
The (almost sure, real) linearity of (2.23) in / guarantees that 

E^{-) = (2.26) 

is indeed an expectation functional (satisfying (E1)-(E3)). Eurthermore, 

RexM)= [ dB^{k)Re{Siik)fik)}- [ dB^{k)S 2 ik)lmfik) (2.27) 

jRd JlRd 


is the sum of two independent normal random variables with mean zero and variances 
||Re 5 i /||2 and || 52 lm/||| (Ito isometry, see below), whose sum equals cj^(/)^ (use (2.24), 
(2.25)). Therefore, ReXcoif) ~ ^/(^(O, cj^(/)^), as desired. This shows the following 
result. 
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Theorem 2.4 Suppose that /i(l) = 0. 


(1) For all f G L‘^{R<^,dk), 

RexM) ~ AA^(0,a^(/)2), (2.28) 

where cr^{f)‘^ is given in (2.18). 

(2) Let and be the functionals (2.19) and (2.26), respectively. Then 

we have, for all f G L^(M‘^,dx), 

En,M) EUf), asN^oo. 

(3) £'aj(-) satisfies (E1)-(E3) in the following sense: 

- Eu}{0) = 1 a.e.{uj) _ 

-for all f G L'^{M.^,dx), E^{f) = E^{-f) a.e.{uj) 

- for all K > 1, Zk & C, fk & L^(R'^, dx), k = 1,..., K, we have 

Ylk,k'=i - fk>) > 0 a.e.{u:) 

Remark. For each / we have Ei^{f) G L^(ll,dP). So there is an 11 j C 11 with 
P(Hj) = 1 s.t. E^{f) G C for all u G H/. (That is, we can choose a representative 
of the function which is well defined and finite on a set of full measure.) Given 
fii ■ ■ ■ 1 fx we thus find ..., Hj^, all of full measure, so that l^uiifj) is well 

defined and finite for all u G , again a set of full measure. The latter sum then 

defines again an element in L^(ll,dP). In this sense, we can form finite (or countably 
infinite) linear combinations of Ei^{fj). The set of lo of full measure on which (E3) 
above holds generally depends on the functions fj. 

Construction of Xcoif )- Let / G L^(M'^,dA:) be complex valued. The Ito stochastic 
integral 

[ dB‘^{k)f{k) 

JWd. 

is a random variable on a probability space L^(n,dP). It is constructed for complex 
random variables as for real ones, see e.g. [15]. Let 

M 

4>{k) = ^2:^(/>jj'(/ci) • • • fij^ikd) (2.29) 

Ai=l 

be a finite linear combination, G C, of indicator functions (fji^{ki) ■ ■ ■ 4>i^{kd), where 
k = {ki,..., kd), kj G M and is the indicator function of the interval = [ak,b^) C 
M. The Ito integral of is defined to be 


M d 

11=1 j=i 
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where B'j, j = are d independent Brownian motion random variables on 

L^(17,dP). Using that the Brownian increments ;= — B‘^(aj) are inde¬ 

pendent for different j, as well as for j fixed and different and that they are normal 
with 

E[AB‘^J=0 and E[{AB‘]’^^f] = b>; - a<;, 
one readily verifies that the following Ito isometry holds: 


E 


dB‘^{k)^{k) 


2l 



(2.30) 


This isometry allows us to define the Ito integral for a general / G (M.'^, dk) as 
follows. Take any sequence (/>n(k) of simple functions (2.29) s.t. (/>„—)■/ in L^(R'^,dk), 
as n —k oo. By (2.30), the sequence of random variables dB‘^(k)(/>n(k) is Cauchy 
in L^(n,dP). It thus converges in the L^(n, dP)-sense to a limit which we take as the 
definition of dB‘^{k)f{k), 


[ dB‘^ik)f{k) 

jRd 


lim 

n^oo 


dB^^{k)(j)n{k) in L^(M‘^, d/c)-sense. 


The Ito isometry then extends to all / G L^(M'^,dA:), 


E 


dB‘^{k)f{k) 


2l 



(2.31) 


2.4 Hilbert space representation 

In the previous section, we have constructed the following expectation functionals of 
the Weyl algebra in the infinite volume limit. 

• describing N discrete modes, see Proposition 2.1, (2.7), 

• {E){f) = limTv^ooE[£' 7 V,i.d(/)] describing the phase averaged functional in the 
continuous mode limit, where is E]\f in which the phases are considered to 
be i.i.d. random, see Proposition 2.2, (2.17), 

• Ei^{f) describing independently distributed random phases in the continuous 
mode limit (with convergence in distribution), see Theorem 2.4 and (2.26). 

Since ReXaj(/) is normal with mean zero and variance cr‘^{f) (see (2.18)), its char¬ 
acteristic function is given by = e“ 2 * o'm(/) _ Therefore, we have from 

Proposition 2.2, (2.17), that 

{E)if) = E[EM)] = f dP{u;)EM)- (2.32) 

Jn 

Relation (2.32) will yield a representation of the state associated to {E){f) as a direct 
integral over the representations of E^^^f). Furthermore, (2.32) shows that taking the 
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expectation E and the continuous mode limit N ^ oo are commuting operations, in 
the sense that the following diagram commutes, 

E 

En,uj -> ^[En,uj] 


N \ I iV (2.33) 

i I 

E 

E^ -^ {E) 


The left down arrow has is a limit in the sense of distributions of random variables. 

Given a state p on a C^-algebra 2t, there exists a unique GNS triple (^,7r, T) 
consisting of a Hilbert space 7i, a representation map tt : 21 —)■ B(Ti) and a normalized 
vector ^ gT-L, such that for all H G 21, 

p(H) = (T,7r(H)T). (2.34) 

Both Etv and (E) define states on the Weyl algebra, determined by pi{W (/)) = E]\f{f) 
and p2{W{f)) = {E) (/), / G L\R<^,dx). 

Gonsider now the family E^. Being an element of L^(H,dP), E^{f) is only well 
defined (represented by a function with finite complex values) for cu G C H for some 
with P(Hj) = 1. The range over which cu varies thus depends on /. Therefore, it 
is not clear that there is any uj G Q for which one can define simultaneously Ei^{f) for 
all / G L^(M'^,dx). However, we can restrict the range of / to a countable subset of 
“test functions” E C L^(M'^,dx). For each f G E, there is an C H, P(Hj) = 1, 
on which Ei^{f) is well defined, that is, for which one can choose a representative of 
the L^(H,dP) function which is finite on Hj. Being a countable intersection of sets of 
measure one, the set 

H(P) = r\f^x> 

has also measure one. Furthermore, for every oj G Q{E) fixed, Ei_j{f) is well defined 
for all / G P. Without loss of generality, we may assume that P is a vector space over 
the countable field Q + iQ. (If an original E is not, then we can consider the set of all 
finite linear combinations of elements of E, with scalars from Q + iQ. This is again a 
countable set and we can take that set for E.) 

We point out that it follows from the construction of the ltd integral that Vcu G H(P) 
and for all / G P, z G C, Xuii^f) = zxuiif), understood as an equality of finite complex 
numbers, for suitable representatives of the L^(P, dP) objects. However, full linearity of 
Xui, using complex scalars, cannot be guaranteed to hold for all uj G fl(E). Nevertheless, 
it follows from the almost everywhere linearity of the ltd integral (see [15]) that for 
every oj G H(P), / Xuiif) is (Q + iQ)-linear on P. (This means that VA, p gQ + iQ, 
V/, p G P, we have Xuii^f + w) = and the last equality holds 

Vw G H(P), for suitable (finite) representatives of the Xaj(')-) 

It follows from the above discussion that there exists an H(P) C H of full measure, 
such that Vw G H(P), is an expectation functional on the Weyl algebra with test 
functions / G P, i.e., the Weyl algebra (over the field C of scalars) generated by all 
W{f),fGE. 
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Theorem 2.5 (GNS representation) 


(1) The GNS representation associated to Efq{-) is given by 

n = L‘^{S^ X ■■■ X S^,dai---daN) 

vrjv(VF(/)) = 

= r^Fock ® 1- 

Here, dcr is the uniform measure on the circle and 1 is the constant function. 

(2) The GNS representation associated to (E) (•) is given by 

n = 'HpCE{L\M.^,dx))^E{L^{R^,dx)) 

T^p{W{f)) = i^Fock(i?/)®T^Fock(r/) 

'I' = ^^Fock ® ^Fock- 

The maps R,T : L‘^{W^,dk) —)■ L^(M'^,dA:) are real-linear, given by 

{Rf)ik) = ^/l + pik)aik)fik) + ^mm (2.35) 

{Tf){k) = ^/l + p{k) P{k)f{k) + y/pa{k)f{k), (2.36) 



(2.37) 

(2.38) 


For fi{2) = 0, we have a = 1 and /? = 0. 

(3) LetV C L^(M^,dx) be a test function subspace with associated Q(V) satisfying 
P(17(D)) = 1. For every oj G FiifD), the GNS representation of E^{-), as a 
functional of the Weyl algebra with test functions in V, is given by 

n = nvQiF{L^{R‘^,dx)) 

TT^iWif)) = iTFock(/)e‘^"^“(^) 

= f^Fock- 


A proof of Theorem 2.5 is obtained by direct verification. In finding the GNS 
representations, we were of course inspired by the Araki-Woods representations, see [1] 
and also [14] for a textbook explanation. 

A representation tt is called regular if a i-G iT{W{af)) is differentiable at a = 0, in 
the strong sense on a dense domain in H. For regular representations, one defines the 
represented Weyl operators 

WM) = Aw{f)) 
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and the represented field operators by 


^M) = -ida\a=o'n:{W{af)) 

and similarly, the creation and annihilation operators by 

alif) = (2.39) 

aM) = 2-V2[$^(/) + i$^(i/)] =(a;(/))*. 

It is not hard to see that the three representations of Theorem 2.5 are all regular, and 
that the creation and annihilation operators are given as follows. 

Proposition 2.6 (Field and creation operators) 

(1) The field and creation operators associated to are 

N 

^wif) = ^*Fock(/) (8) ]1 l 2 - Ijr (g) ^ |cos6'jRe/(%) + sin9jlmf{kj)'j 

i=i 

N 

i=i 

(2) The field operators associated to (E) (•) are 

^p{f) = $Fock(R/)®l.F + ]l.F®fRFock(r/). 

Since R, T are only real linear, the creation operators have a somewhat cum¬ 
bersome expression. For fi{2) = 0 it reduces to 

o*p{f) = apockiV'^E Pf) (8) 1 + 11.F ® apockiy/pf)- 

(3) The field and creation operators associated to E^, for all u G Ft.{V) and all 
f GV, are 


^M) = ^Fock{f) + RexM) 

«w(/) = «Fock(/) + :(^ Xw(/)- 


Let fi, ..., fn G V. For every lo G fI(P), we have (strongly on a dense domain of 
Hilbert space, for instance on the finite-particle subspace of Fock space [4]) 




i-iy 

{-iy 


QT. 


dai ■ ■ ■ don 

I 

dai ■ ■ ■ don 


WUaifl)---WUotnfn) 
e'^fLaj(ai/l + • • • + Onfri), 


where P is a real, deterministic phase depending on the aj and the fj (originating from 
the Weyl canonical commutation relations (2.1)). 
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Denote by (•)^ the average in the state defined by the vectors 'h in the GNS repre¬ 
sentations (in the appropriate Hilbert space). Relation (2.32) then implies that 




dai 

■ ■ ■ dan “1“ 

■\n 

d^ 1 

dai 

• • • dan "1= 

■xn 

d^ 1 

^ dai 

• • • dan "1= 

Afi)-- 

• ^p{fn))q, ■ 


.=«n=0 

.=0,1=0 

.= a „=0 


e'^ (E) (ai/i H-h a„/n) 

e'^(Wp(ai/i H- ^anfn})q, 

{Wpia,fl)---Wp{anfn))^ 


(2.40) 


Relation (2.40) can also be expressed in terms of creation and annihilation operators 
as follows. Using (2.39) and forming linear combinations of (2.40) we see that 


E 


^a;(/l) • • • ^Mn-l)a*{fn))^] = (^p(/l) ' ' ' ^p(/n-l 



where stands for either a or a*. Then we can continue the procedure to see that 
in the relation (2.40), the product of held operators can be replaced by a product of 
arbitrary creation and annihilation operators. 

This allows us to show that n point functions can be expressed in terms of two 
point functions only. For fi,..., fp and gi,... gq p + q test functions in V, we dehne 
the block matrix Q G Mp+g(C) by 


Q 


A 

C B 


(2.41) 


where A G Mp{C), B G Mg(C), C G Mg^p(C) are dehned by 

Aij = '[l{2){fi\pfj), Bij = p{2){gi\pgj), Cij = {gi\pfj) 
and where is the transpose of C. 


Proposition 2.7 (Quasifreeness) Let fi,..., fp and gi,... pg be test functions in V. 
Then 


IE[(a()(/i) • • • al{fp)au{gi) ■ ■ ■ aui{gq))-i,] 

Jo if p + q is odd 

\ XjttSGJj Q7r(l)7r(2) Q7r(3)'7r(4) • • • Q'K{n—l)'n-{n) ^ P q is CVCn. 


(2.42) 


Here, is the set of all permutations of (1,... ,n) such that 7r(l) < 7r(3) < • • • < 
7r(n — 1). In particular, if p{2) = 0, 


■ ■ ■ al,{fp)au,{gi) ■ ■ ■ 



0 ifpf^q 

EaeSp ]Tj=i{9a{j) \pfj) ifp = q- 
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In particular, we have 


maUfi)h] 

maUfi)a*M2)h] 

IE[(aa;(5fi)a^(5'2))4'] 

E[(a^(/i)aa;(5'i))4'] 

E[(a(^(5'i)a^(/i))4'] 


E[(a,(/i))vi,] =0, 
Ai(2)(/i,/5/2) , 
m ( 2 ) {gi,pg2), 

{gi,pfi), 

{gi,{p + l)fi) . 


Note that the functional (2.42) is not gauge invariant unless p{2) = 0. I.e., it has 
nonzero average of products with unequal numbers of creation and annihilation oper¬ 
ators if ju(2) / 0. 


2.5 Dynamics 

2.5.1 Reservoir dynamics 

The dynamics on the Weyl algebra is given by a Bogoliubov transformation on the 
functions / G L^(M'^,d/c), 

/ ^ 

where s = e{k) is a real function of /c G As an example, for photons, e{k) = |A:|. 
The dynamics of the three infinite-volume expectation functionals is given as follows. 
For N discrete modes, (2.7), 

= i^Fock(/) (2.43) 

where we note that Epocki^'^^ f) = -EFock(/) for all t and all /. The dynamics has thus 
the effect of shifting the phases associated to the coherent states. 

The phase-averaged expectation functional (2.17) evolves according to 

(E) (e'*^/) = Ep,,p{f) (2.44) 

where 

= ^^p(A:)(^|/(A:)|^-bRe{e^‘*^/i(2)/(/c)^})d/c. (2.45) 

It follows from the Riemann-Lebesgue Lemma that 

hm u^(e‘*7)2 = W^fWl (2.46) 

r —¥ cxD 

To see how to derive (2.46) in the situation where e = £(|A:|) is an invertible function 
of |A:|, with inverse \k\ = v{e), we use spherical coordinates in W^, 

p poo 

/ e^^^p{k)f{kf = / e2‘*^F(y)dy, (2.47) 

Jo 

where E{y) = v{yY~^v'{y) f^d-i J0) p{v{y),6)f{v{y),9)‘^d9. Then the ordinary Rie- 
mann Lebesgue Lemma asserts (2.46) provided E G L^(M+,dy). 
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Therefore, the infinite volume state converges, for large times, to the state deter¬ 
mined by the uniform phase distribution (d/i(0) = ^). Note that any phase distribu¬ 
tion jj, satisfying ^(2) = 0 determines also that same state. 

The random phase expectation functional (2.22) satishes 

= ii;Fock(/) (2.48) 


Due to (2.46), 


A/Li(0, (T^(e^*^/)^) )■ A/’aj(0, IlyT^/lll), t —)■ oo. 

Therefore, the random infinite-volume state converges to the the random infinite volume 
state with uniform phase distribution in the limit of large times. This convergence in 
the sense of distributions of random variables. 

We have shown the following result. 

Proposition 2.8 (Phase uniformization under reservoir dynamics) Let be a 

phase distribution satisfying = 0. Given any f G L^(M'^,dx), we have, as t ^ oo, 

(E) (e^‘V) ^ (i^)unif (/) (2.49) 

A £;.;,unif(/). (2.50) 

The convergence in (2.50) is in distribution of random variables. Here, (£^)^j^j£ (•) 
and £'a;,unif(') the expeetation functionals (2.17) and (2.26) in which the phase 
distribution is uniform, dfj,{9) = 

Remark. All measures d/i(0) satisfying j(i(l) = 0 and having the same value of ju(2) 
give the same expectation functional {E) and E^ (see Proposition 2.2 and Theorem 
2.4). 

2.5.2 Coupling to an open quantum system 

We consider an A^-dimensional quantum system in contact with the reservoir of coherent 
states in which the phases are uniformly randomly distributed. The Hilbert space of 
pure states of the system is , that of the reservoir is the GNS space given in point 
(3) of Theorem 2.5. The system dynamics is generated by a self-adjoint Hamiltonian 
with energy levels ei,..., ev, 


Hs = diag(ei,... ,eAr). 

The state of the reservoir is invariant under its own dynamics by Proposition 2.8. The 
dynamics is implemented as 

7rUW{e^^^f)) = lHFock(e‘'"/)e'^'^>^“(4) = e'*^^7r^(W(/))e-'*^^ 
where the reservoir Hamiltonian is 


Hr = dr(e). 
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The uncoupled dynamics is therefore given by the Hamiltonian 

Ho = Hs 0 + ts ® Hr. 

To define a coupled dynamics between the system and the reservoir, one proceeds as 
follows. The free dynamics is given by the group of ^automorphisms Og on the algebra 
of observables 2t = H(C^) ( 8 * W (where W is the Weyl algebra), defined by 

(8) W{f)) = <8) 

Then one defines a coupled dynamics by specifying an interaction operator H G 21 and 
using the Dyson series 

rt j-tn — l 

a\A) = al{A)+ Y, dti • • • / dt„ [«*" (H), [• • • a* (H)] • • • ]]. (2.51) 

n>l -^0 do 

The series converges in the topology of 21 and defines the interacting dynamics a*, again 
a grop of ^automorphisms on 21. Applying the representation map Tr^; (more precisely, 
Is ® to (2.51), we obtain 

T^Lj{oi\A)) = To{tt^{A)) (2.52) 

pi rin 1 

+ ^ / dti--- / dfn [ro’*(7r^(H)),[---[roi(7ra;(H)),ro*(7r^(A))] •••]], 

n>l -^0 -^0 

where 

r*(.) = e^*-^°(-)e-^*-^°. 

The right side of (2.52) defines a group of ^automorphisms on the represented algebra 
of observables which is generated by the self-adjoint operator 

H = Ho + ti^{V), 

acting on C'^ ( 8 > T'(T^(M^, dx)). From physical considerations, one would like to take 
V = G(S>^(g), where G is selfadjoint and ^(g) is a field operator. Of course, this V does 
not belong to 21 and the above construction cannot be carried out. Nevertheless, one 
can “regularize” the interaction by introducing V^, depending on a small parameter g, 
such that G 21 and in any regular representation tt of the algebra 21, 7r(F?) —)• G( 8 i<h 7 r, 
as ry — 7 - 0 (strongly on a dense domain). One can then, for g > 0, carry out the above 
construction and finally remove g once placed in a representation. Such a procedure 
is decribed in [5] - and other approaches are possible. The dynamics of the coupled 
system is thus generated by the Hamiltonian 

H = Hq + G ® d>t^( 5 ') = Hq + G ® (d>Fock(< 7 ) + ReXa;(fl')), (2.53) 

acting on ® F{Lp‘{W^, dx)). 

We consider an energy eonserving (non-demolition) interaction [10,16] between the 
system and the reservoir, which consists in taking an operator G that commutes with 
Hs, 

G = diag{gi,...,gN)- 
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Such models are used to investigate “phase decoherence” of the small system. 

The initial system-reservoir state is disentangled, given by a density matrix 

Pq = PS® |^)(^|, 

acting on the Hilbert space C'^ (8>dx)). Here, ps is an arbitrary intial system 
density matrix and the reservoir is in the state H, which represents the infinitely ex¬ 
tended continuous mode coherent state with uniformly distributed phases. The state 
of the coupled system at time t is given by 

P(t) 


Taking the partial trace over the reservoir Hilbert space yields the reduced system 
density matrix, 

Ps{t) = TrRP(t). 


We denote its matrix elements in the energy eigenbasis 
by 


Pk,l{^) = {pk,ps{t)pi) = P{t)\Pl){Pk\- 


(with H^ipj 


= em) 


(2.54) 


As [Hs,G] = 0 the populations (diagonal matrix elements) are time-independent. The 
off-diagonal ones exhibit time decay (“phase decoherence”). For the energy conserving 
model at hand, the matrix elements (2.54) can be evaluated exactly. The calculation 
yields (see Appendix D of [13]) 


xpk,i{0)- 


(2.55) 

(2.56) 


The contribution on the right side of (2.55) is given by the free dynamics and by a 
random “renormalization” of the system energy due to the interaction with the coherent 
bath (coming from the term G(8>ReXa;(5)]lR in the Hamiltonian (2.53)). The two factors 
(2.56) are the same as if the system was coupled to a free bose gas in equilibrium at zero 
temperature. Therefore, coherent states character of the reservoir is encoded entirely in 
the part e~^^Gk-gi)^exujig) ^ expectation of this oscillating factor is the characteristic 
function of the random variable ReX(.j(5') MU0,\\Vp9\\1), 


E 


g-ii(9fc-9i)ReXuj (g) 


= g-4(Sfc-3i)^ll\/P9ll2, 


This shows that the averaged (reduced system) density matrix E[ps(i)] acquires Gaus¬ 
sian time-decay of off-diagonals at all times, due to the coupling with the coherent 
reservoir, namely 


with 


|lE[pfc,z(t)]| = e-'^^3k-gi?\\^pg\\l^-\{gu-gi?nt) \pj^i{0)\, (2.57) 



1 — cos(ef) 


sin^(et/2) 

17) l9 9 


(2.58) 
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For small times, T{t) ~ ^^^115112 is quadratic in time, but for large t, its behaviour as a 
function of t depends on the infrared behaviour of the form factor g. For instance, in 
d = 3 dimensions and for £{k) = \k\, 

T{t) = 2 r\k\Mk\ [ hm [ dS|5(r,S)|2, 

Jo Js^ Fr 2 r^0+ Js2 

assuming that the latter limit exists and is non-vanishing, meaning that | 5 '(r, S)| ~ 
for small r. Note also that for p > —1/2, we have limt^oo {d, cos{st)/£‘^g) = 0 by the 
Riemann-Lebesgue lemma, so that limt^oor(t) = Us'/eHl- For this infra-red behaviour 
of the form factor, the coupling to the (zero temperature) reservoir does not induce 
(complete) decoherence, but the coupling to the coherent reservoir does. 


3 Proofs 


3.1 Proof of Proposition 2.1 

We calculate the limit L —)■ oo of 

EUf) = 




W(/)$) 


(3.1) 


where W (/) = = e fka.k+fk<^k Weyl operator in momentum space. The 

coherent state (2.3) has the form 

$ = 1-09)0, with = 

and where the momenta satisfy limi^oo (T) = kj and aj{L) is given in (2.6). Using 
the canonical commutation relations W{f)W{g) = e~ + g) yields 

E^f) = {n,W{-g)W{f)WCg)n) = 

Combining (2.6) with the formula (2.2) gives 


N N 

'^aj{L)fkUL) = 

i=i i=i 




f{x)dx, 


which converges to \/7h® /r'* ® ^^^'^/(x)dx in the limit L —)■ oo. 


□ 


3.2 Proof of Proposition 2.3 

The mechanism of the proof of Proposition 2.3 is the following. Since the defined 
in (2.20), are independent with mean zero and variance 

EK|| = {2R)'‘p{ki)[\hki)f + Re{9(2)/(t:j)n), 
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(3.2) 


where Jl{2) is given in (2.18), a version of the central limit theorem says that 


.A V„(0,1), 

jv-iEjEd. K)^n(] 


as —)■ oo. 


Here, A/'aj(0,1) is the standard normal with mean zero and variance one. We point out 
that the are not identically distributed (only the 6j are). We give a proof of (3.2) 
here below. The denominator in (3.2) has the limit cr^{f) given in (2.18), as —)• oo. 
Hence (3.2) implies 


N ^ crA‘(/)A/L(0,1) = A/L(0,o-^(/)^), asA^-^oo. 


We now prove (3.2). Define 


s%:= p{kj)[\f{kjf+Re{mmf})- 

The Central Limit Theorem (see Theorem 27.3 in [3]) says that 


1 


— ^j(^) A^w(0,1), as A^ —)■ oo, 

C^r ^ 


sn 

provided that the following Lyapounov condition holds, 

E[ie,(v)|2+'] 


lim 

IKJ _i ("V"! 


N^OO ‘ ~~ Q.^ 

je{i,...,NY ^ 


2+(5 


= 0 


(3.3) 


(3.4) 


for some <5 > 0. We have 

IE[|,ej(w)P+'5] < (2i?)^(2+5)/2 (2 /j(A:j))(2+'5)/2|/(A;^.)|2+5 ^ 

as well as 

< (2i?)'^(2+^)/2| Y Mkj) ~ Ar'=*(2+5)/2^ 

Therefore (3.4) holds since the sum is of the order . □ 


3.3 Proof of Proposition 2.7 

Consider the complex valued random variable 


{a*Mi) ■ ■ ■ a*Mp)a^ 


iai) 

1 




.+qXMl) ■ ■ ■ XMp)Xaj{9l) ■ ■ ■ xU9q) 


(3.5) 
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For functions /i, • • • , /p, ffi, •'' > complex numbers zi, - ■ ■ ,Zp,wi 

define 

$(z, re) = E exp (E ^jXuiifj^ “1“ ^ '^kXuijdk)} 

j k 

The derivative at the origin yields the sought for expectation value, 


, ••• ,Wq, 


= E xMi) ■ ■ ■ xMp)Xio{gi) ■ ■ ■ xUdq 


We have 


'^ZjxMj) + '^WkXu,{gk) = XioC^Zjfj)+XujC^Wkgk) 


(3.6) 


(3.7) 


(3.8) 


= Xuj,. 


dE + '^Wkgk] +iXuj,2 Wkgkj 


where Xui,ii') Xa;, 2 (') are the real and imaginary parts of Xu){')- 

Consider the characteristic function of (Xaj,i(/))X i.j, 2 ( 5 '))- For any real ri,r 2 , the 
random variable riXi.j,i(/) + T 2 Xui, 2 {g) is normal with zero expectation value and 

(3.9) 


-\^^QK'riXu,,l{f)+T2Xuj,2 

= exp [rfyaiXioxif) + r|varxa;,2(5') + 2rir2Cov(x<.,;,i(/),X(.^,2(5'))] 


Note that ti,T 2 can be taken here as arbitrary complex numbers, since the distribution 
of Xuiji'), j = 1,2, is Gaussian. 


Lemma 3.1 With the notations above, 


^{z, w) = exp < 2 E ZjWk{gk\pfj) 

( j,k 

+ /1(2) ^ ZjZf{fj/\pfj) + fi{2) ^ wkWk'{gk'\pgk)\ ■ (3.10) 

j,j' k,k' J 

Proof: According to (3.6) and (3.8), is given by (3.9) with ti = — i, r 2 = 1 
and / = Y^jZjfj + Y^k^kgk, g = JljZjfj - Set a = Yfj=iZjfj and b = 

^k=i^9k- Using (2.23), (2.24) and (2.25), we compute 

varxi.j,i(a+fc)-varXaj, 2 (a-&) = J pik) ^2{ab + ab) + A(2)(a^ + b^) + p{‘2){a? + 62)| dk, 

and 

cov(xa;,i(a+^),Xa;, 2 (a-&)) = ^ J |2(a6 - db) + A(2)(a^ - b^) - p{2){a^ - 62)| dk, 

from which (3.10) follows. □ 
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We use the notation t = {z,w) G x and write accordingly $(t) = 

^{z,w). Lemma 3.1 shows that <h(t) = where q{t) = n}^j^kQjk is a 

quadratic form, with corresponding matrix Q G Mp^g{C) defined by (2.41). Equation 
(3.7) takes the form 


E 


xMi)--- xMp)xMi) ■ ■ ■ Xcoigq) 


= d: 




E 

e>o 


ei 


( tjtkQjk) 

j,fc=i 


(3.11) 


The right hand side is not zero only if re is even and only the term i = nj^ in the series 
does not vanish, so that we have 


E 


Xu){fl) ■ ■ ■ Xuj{fp)Xoji9l) ■ ■ ■ Xioidq] 


- ^*1,. 




(re/2)! 


^ ^ tjtkQjk) 
j,fc=l 


n/2 


(3.12) 


The power re/2 of the double sum over j,k in the right hand side of (3.12) is an re¬ 
fold sum over values tg, s = l,...,re and each summand has a factor of the form 
tjitki • • • only terms in that multiple sum that have a nonzero derivative 

are those in which (ji, ki,... ,jn/ 2 i ^n/ 2 ) are a permutation of (1,2,..., re). Therefore, 


E 


w(/i) • • • xMp)xU9i) ■ ■ ■ Xw{gq) 


1 

(n/2)! 


^ ^ Q'k{1)'k{2) ' ' ' Q'K{n—l)'K{n)- 

TT&&„ 


(3.13) 


Each term in the sum of (3.13) has the same value upon permuting the re/2 different 
factors Q, so we arrive at (2.42). This concludes the proof of Proposition 2.7. □ 


3.4 Simultaneous continuous mode and infinite volume limits 

We consider the one-dimensional case, d = 1. The discussion for general d is analogous. 
Suppose we want an infinite volume coherent state with modes distributed in an interval 
[a, b] C M. The finite volume state is 

The discrete modes in [a, b] are kj = a + j^{b — a) and the sum in the exponent is 


L/n 

^ ^ ^kj 0-kj ~ ^kj O'kj ■ 
1=0 


This gives the expectation functional (see (2.7)) 


Ehig) = EYocv{g)e 


iv^Re E/ii 9kj 


£^Fock(fl')e 


E/ii gikj) 


(3.14) 


The total particle density is p = ^ = Y^^^Pkj- Let p{k) be the density 

distribution, i.e., fj p(k)dk is the density of particles having moment in the interval 
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I C [a,b] C M. Then p{kj) = pkj/Akj = p^.L/tt and hence p = p{kj)/^kj, so 

p = p{k)dk in the limit L —)■ oo. The link between ak^ and p{kj) is 

oikj = “^Lpkj = yj2Lp{kj) \/A^ = yj2Txp{kj) e‘®''> , 
since Afcj = vr/L. We use this in (3.14), 

Edg) = i^Fock(l/) 

The phase behaves like 2^J L/tt Re p{k)e~^^^^^g{k)(\k and diverges like \/L. 


3.5 Rarefied continuous mode limit 


We consider d = 1. General d are treated in the same way. The origin of the non¬ 
existence of the continuous mode limit (2.9), (2.10) is the following. The discrete 
density is related to the continuous density distribution p{kj) by pkj oc -^p{kj), 
where 1/A is the discretization mesh size. The phase of the expectation functional is 
a sum over the square root of the density, Ylf=i y/TEJoi^j) — Si=i VPi^j)gi^j)^ 

giving a prefactor which is too weak for convergence of the Riemann sum. 

Instead of populating all modes, we may consider a rarehed situation, where, in 
a given interval [a,b], only a fraction of modes are chosen. The modes in the entire 
interval are kj = a + j^{b — a), j = 0,..., L/tt. Fix an integer s and choose 

f (Ufc ■ ior j = is, £ = 0,-E 
^ \ 0 for all other j 

The sum in the phase of (3.14) becomes 


L/'n 


^ X/ 

^ jr = l 


^ L/sn 

j-jr ^ ^ ^a+ts'n{b—a)/L fl'(® T is 7 r(^b <x)/L). 

vT 


The mesh size in the Riemann sum over i is A/c = ^{b — a), and we see that the right 
side converges, and has a nonzero limit, only if {Aky/L)~^ converges to a nonzero limit. 
This forces s oc y/L. Taking s = for some u > 0, we obtain 


lim 

L—^-OO Jj 



a{k)g{k)dk. 


The infinite volume expectation functional becomes 

lim Eiig) =F^Fock(d) . 

L^oo 


We examine the particle number and density, as L —)• oo. The number of particles is 


Lj-K L/s-k . 

= l«a+fo7r(6-a)/Lp ~ / \a{k)\‘^dk. 

j=o e=i 
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(Note, the mesh size in the last Riemann sum is AA: = ^{b—a) = {ay/L)~^.) Therefore, 
the number of particles grows like a times the square root of the volume. Dividing by 
the volume, we obtain that the total particle density decays to zero like We 

have thus an infinite volume limit state containing infinitely many particles, but having 
zero particle density - a rarefied state. 
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